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We present an exploratory calculation of the i" = 2 tt — tt scattering amplitude at 
threshold using Wilson fermions in the quenched approximation, including all the required 
contractions. We find good agreement with the predictions of chiral perturbation theory 
even for pions of mass 560-700 MeV. Within the 10% errors, we do not see the onset of 
the bad chiral behavior expected for Wilson fermions. We also derive rigorous inequalities 
that apply to 2-particle correlators and as a consequence show that the interaction in the 
antisymmetric state of two pions has to be attractive. 
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1. Introduction 

In this paper we calculate the 1 = 2 tttt scattering amplitude at threshold using Wilson 
fermions. The theoretical foundations of the calculation have been established in a series 
of papers by Liischer where he shows how the finite volume dependence of the two particle 
energy levels in a sufficiently large cubic box is related to the scattering amplitude [1] [2]. 
In the context of lattice field theories, the leading term in an infinite volume expansion 
was first given in Ref. [3]. The subtleties of the lattice calculations have been described in 
detail in Ref. [4], where results of a calculation using staggered fermions is presented. We 
follow closely the notation of this reference. 

The calculation of the 1 = 2 tttt scattering amplitude at threshold has a number of 
simplifying features. In general, scattering amplitudes are complex, and are only related 
indirectly to the finite volume energy shift. At threshold, however, the amplitude is real, 
and is directly related to the energy shift. Also, the signal for pions is much better than 
for other, heavier, mesons, e.g. rhos. Finally, one needs to calculate only quark and gluon 
exchanges for the 1 = 2 channel, whereas, in general, there are also annihilation diagrams. 
The latter are both more difficult to calculate numerically, and are affected more strongly 
by the use of the quenched approximation. 

A major motivation for this work is to test the chiral behavior of the scattering 
amplitude derived long ago by Weinberg [5] using PCAC and current algebra. We use 
Wilson's formulation of lattice fermions, for which there are lattice artifacts arising from 
the explicit breaking of chiral symmetry. We can quantify these corrections by comparing 
lattice results against the PCAC prediction and against results obtained using staggered 
fermions on the same set of lattices [4]. 

Our calculation is exploratory in at least three different ways. First, we use the 
quenched approximation. Second, we only use two different quark masses, neither very 
light. Thus we can only make rough extrapolations to the chiral limit. And, finally, we 
have only used one lattice size (16 3 X 40). This means we must assume the finite volume 
dependence predicted by Refs. [1] and [3], and cannot check the predictions. 

In the earlier calculation, Guagnelli, Marinari and Parisi [6] have done a partial cal- 
culation of the scattering amplitude using Wilson (and staggered) fermions. They did not 
include all the contractions which contribute in the 1 = 2 channel, and thus one cannot 
extract the scattering amplitude from their results. 

The article is organized as follows. A brief theoretical overview is given in Section 
2 and the methodology and details of the lattices are given in Section 3. The results 
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are presented in Section 4, and in Section 5 we present a model which explains some 
features of the data. Section 6 gives some conclusions. In the Appendix we derive rigorous 
inequalities regarding correlators of the type used in this study, in particular we show that 
the interaction in the anti-symmetric tt — tt channel has to be attractive. 

2. Theoretical background 

Liischer has derived the relationship between pion scattering amplitudes and the en- 
ergies of two pion state in finite volume [1][2]. The derivation is valid as long as the box 
(which we take to be cubic) is large enough that its length L exceeds twice the range of 
interaction. In general, the relationship is complicated, involving scattering amplitudes 
over a range of energies and in many partial waves. The relationship simplifies, however, 
if one expands the energies in powers of 1/L and keeps only the first few terms. We con- 
sider only the lightest two pion state whose energy we denote by E. For infinite volume 
E = 2m K , but the energy is shifted by interactions as L is reduced 

6E = E — 2m n = (l — d -j^- + c 2 ( -^-X) + 0(L" 6 ) , 

L 3 \ 47rL \AttL J J (2.1) 

ci = -2.837297 , c 2 = 6.375183 . 

To the order shown, the energy shift depends only on T, which is the (non-relativistically 
normalized) scattering amplitude at threshold. 

We use the non-relativistically normalized amplitude in Eq. (2.1) since this simplifies 
its physical interpretation, as explained in Ref. [4]. T is related to the relativistically 
normalized scattering amplitude by T R = — (4m^.)T, and to the S-wave scattering length 
by T = — 47rao/m 7r . 

We extract T from our numerical data using Eq. (2.1). A priori we do not know the 
size of the 0(L~ 6 ) terms which we are dropping. Our numerical results suggest, however, 
that the truncation error is small. 

Equation (2.1) holds separately for 1 = and 2 two pion states (Bose symmetry 
forbids / = 1 at threshold). We have done the calculation, however, only for the 1 = 2 
channel. To understand why, consider the four types of contraction that contribute to a 
calculation of the two pion energy, shown in Fig. 1. With present computer resources 
we can calculate only the first two types, which we label the Direct (D) and Crossed (C) 
diagrams respectively. These are not sufficient to calculate the 1 = amplitude which gets 
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contributions from all four diagrams. For 7 = 2, however, quark-antiquark annihilation 
is not possible, and the required result is given by the combination D — C . Here we are 
adopting the convention of showing explicitly the minus sign from Fermi statistics, so that 
the lines in Fig. 1 represent c-number quark propagators. 

We can also calculate the combination D + C . This does not project onto a definite 
isospin, but does select a definite representation in a theory with Nf > 4, Nf being the 
number of flavors. D + C picks out the ~qqqq representation having no traced indices, and 
is antisymmetric under the interchange of either quarks or antiquarks. We refer to this 
representation as the A (for Nf = 4 it is the 20), and to the corresponding scattering 
amplitude as T(A). For arbitrary Nf > 2, the combination D — C, which we call the 
S, projects onto the representation with no traced indices and is symmetric under quark 
exchange. The generalization of the 1 = 2 amplitude T2 is thus T(S). In the quenched 
approximation the amplitudes T(S) and T(A) are independent of Nf, because the Wick 
contractions are always the same. In particular, T2 = T(S). 

An important test of any calculation of pion scattering amplitudes is that they satisfy 
the constraints of chiral symmetry. In particular, the threshold amplitudes are determined, 
in the chiral limit, in terms of f K (which is 93 MeV in our normalization) [5] [4] 



These results should apply in the quenched approximation, as discussed in Ref. [4]. Since 
T(S) > 0, two pions in an S representation are repelled in the chiral limit, while in the A 
representation there is an attraction of equal strength. This equality can be understood as 
follows. The diagrams of Fig. 1 serve dual purpose. In addition to showing contractions 
contributing to E, they can also represent the contributions to a direct calculation of pion 
scattering amplitudes. We refer to Fig. la as the gluon exchange amplitude, T 9 , and to 
Fig. lb as the quark exchange amplitude T q . Following the standard usage for amplitudes, 
we include the Fermi-statistics sign in T g , i.e. we use the opposite convention to that for 
C . Thus we find that T(S) = T g + T q and T(A) = T g — T q . Now, it is possible to show 
that T g vanishes in the chiral limit [4], so that T(S) = T q = — T(A). 

Testing the relations Eq. (2.2) is particularly important for Wilson fermions, which we 
use here. This is because Wilson fermions explicitly break chiral symmetry, the violation 



4/£T(S) 
4/£T(A) 



1 + 0{m\ mm*-) , 
— 1 + 0(m\ In m n ) 
-4/£T(S) + 0(ml 



(2.2) 
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only vanishing in the continuum limit. Thus we expect T to contain a lattice artifact 
proportional to Aa/m^., where a is the lattice spacing and A is some non-perturbative 
scale. This artifact will dominate over the constant term in the chiral limit [7]. This 
is in contrast to staggered fermions where such artifacts are forbidden by the residual 
chiral symmetry [4]. With Wilson fermions one can isolate the physical result by doing 
the calculation at a number of values of the quark mass provided /^T is a well behaved 
function of m\. This remains to be checked. We find that, for 560 MeV < m w < 700 MeV 
the chiral symmetry breaking effects are smaller than the statistical errors (~ 10%). 

3. Calculational details 

The energy of two pions in a finite box is obtained from the Euclidean correlator 

CUt) = (^20 1 (x 1 ,t)^20 2 (x 2 ,t) S 3 (x 3 ,t = 0) S 4 (x 4 ,t = 0)) • (3.1) 

Xl x 2 

The sources Si create the pions at t = 0, and the operators Oi (which we also call the 
"sinks") destroy them at time t. The representation of the two pion state is determined by 
the flavor of the sources and sinks. For example, we can select the 1 = 2 (or equivalently 
S) representation if both sources have the flavor of a tt + , and both 0\ and 2 have the 
flavor of a tt~ . 

At large \t\ the correlator will fall as 

C^(t) = Z^exp(-E\t\) + ... , (3.2) 

where E is the energy of the lightest two pion state. The ellipsis indicates contributions 
from excited states that are suppressed exponentially. This is similar to the behavior of 
the two point function used to calculate 

C„(t) = (J20(x 1} t) S(x4,t = 0)) = Z 7T exp(-m 7T \t\) + ... . (3.3) 

Xl 

We take all quarks to be degenerate, so the flavor of the pion source S in this equation 
is unimportant; all that matters is that O has the conjugate flavor to S. It is useful in 
practice to combine Eqs. (3.2) and (3.3) 

m = ^Jtf = ex P("<® l*D + ■ ■ ■ > ( 3 - 4 ) 
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and directly extract the energy shift SE. 

The contractions which can contribute to C KK are shown in Fig. 1. The combination 
which is needed depends on the flavor of the two pion state. It is easy to see that only 
the Direct and Crossed contractions contribute for two pions in S or A representations 
[4]. We adopt the notation that D(t) is a ratio as in Eq. (3.4), with the numerator being 
the Direct contraction. Similarly, C(t) is the ratio with the Crossed contraction in the 
numerator. Then, as discussed above, and explained in more detail in Ref. [4], we can 
extract the energy shifts for the S and A representations using 

D(t) + C(t) = Z A e ~ 6EA ^ , (3.5) 

D(t)-C(t) = Z s e- SEs ^ . (3.6) 

The amplitudes Zs and Z a are shorthand for the ratios Z^^jZ^. 

We are free to choose the form of the sources and sinks, as long as they both couple to 
a two pion state of the required flavor. This choice does not affect the value of the energy 
shift, but it does alter the signal to noise ratio. In order to improve this ratio, we should 
use sources with a large overlap with the lightest two pion state. In fact, as indicated 
in Eq. (3.1), we use the product of two independent single pion operators, and make no 
attempt to account for the correlations caused by the interactions between the pions. For 
the single pion sources (the Si in Eq. (3.1)), we use wall and Wuppertal quark sources, 
which we have shown to be reasonably effective in producing single particle correlators [8] 
[9]. In addition, we use both pseudoscalar (V = ^75^) and axial vector (A4 = ^7475^) 
operators for each of the sources. There are thus four sources in all, which we label Wp 
(wall with pseudoscalar), W a (wall with axial vector), Sp (Wuppertal with pseudoscalar), 
and Sa (Wuppertal with axial vector). For the sinks (the Oi in Eq. (3.1)) we use local 
operators, with Dirac structure either P or A. Of the various possible combinations of 
sources and sinks we consider only those in which both sources are of the same type, and 
both sinks have the same Dirac structure as the sources. Thus we can label the ratios lZ{t) 
according to the choice of source, i.e. as Wa, Wp } etc. Finally, we always define the ratio 
lZ{t) with the same sources and operators in both the numerator and denominator. This 
is of no consequence for the energy shifts, but does affect the amplitudes Za,s, which we 
discuss further in Section 5. 

We use 35 pure gauge configurations of size 16 3 X 40 generated at f3 = 6.0. We use 
Wilson fermions at two different quark masses, n = 0.154 and 0.155. The corresponding 



5 



pions have masses of about 700 and 560 MeV respectively. The quark propagators are 
calculated on lattices doubled in time (i.e. of size 16 3 X 80), with periodic boundary 
conditions in all directions. We have previously used these lattices and propagators to 
study the spectrum and matrix elements [8] [9], and we list the relevant results in Table 1. 
The value of f K given in Table 1 is different from that quoted in [8] for two reasons. First, 
we use the normalization such that the experimental value is f K = 93 MeV , and second we 
now use the mean field improved value 0.77 for the axial current renormalization constant 
[10], rather than the previous estimate 0.86. The statistical errors in individual data points 
shown in Figures 2-6 are calculated using the single elimination jack-knife procedure. The 
"forward" and "backward" propagators on the underlying 16 3 X 40 lattices give us two 
results for R(t) on each lattice. Since these are correlated, we average them and treat 
them as a single result. 

4. Results 

To display our results, we use the quantity 

8E eS (t) = \n[K(t)/K(t + 1)] . (4.1) 

This "effective energy shift" should reach a plateau of height SE when t is large enough that 
the lightest state dominates. Figures 2 and 3 show 6E e fi(t) for the pseudoscalar operator 
V at n = 0.154, using Wuppertal and Wall sources respectively. The results using the 
axial operator A4 are of poorer quality, and are not shown. There is a clear signal of a 
non- vanishing energy shift. We extract SE by fitting lZ{t) to a single exponential, selecting 
the range of time-slices separately for each channel based on the extent of the plateau in 
the effective energy shift. The solid lines in the figures indicate the fit value over the range 
of the fit, while the dashed lines show the 1 a jack-knife errors. 

It is apparent from Figs. 2 and 3 that the signal for Wall sources has smaller statistical 
errors than that for the Wuppertal sources. We can partly understand this as follows. The 
Wall source produces pions with p = 0, while the Wuppertal source couples to pions having 
all possible momenta. Consequently the Wuppertal source correlators have an additional 
unwanted contribution at small t from an excited state of two pions of equal and opposite 
momenta, each of magnitude p = 2tt / L. For large volumes this state approaches the lightest 
state consisting of two pions both having p = 0, and provides the largest contamination 
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to lZ{t) rather than the states made up of radially excited pions. It may therefore become 
necessary to use wall sources for calculations on larger lattices. 

In Figures 4-7 we show the ratios lZ{t) at n = 0.154, for both V and A4 operators. 
These plots show that there is a contamination from "wrap-around" effects starting at t ~ 
30. One of the two pions can propagate N t — t = 80 — t time-steps backwards, because of the 
periodic boundary conditions. This results in a contribution which is independent of t, but 
suppressed by roughly exp(— m K Nt)/ exp(—2m 7T t) compared to the forward propagation of 
the two pion state. In practice we always fit to time ranges satisfying t max < 26, for which 
we can ignore this contamination. 

The results of our fits, together with the time ranges used, are given in Table 2. 
For both Sp and Wp correlators we fit using the full covariance matrix over the range 
of the plateau. We are unable to do this for the Sa and W a channels, because some of 
the jackknife samples are too noisy. Our results for these channels are obtained keeping 
only the diagonal elements of the covariance matrix, i.e. neglecting correlations in lZ{t) 
between different time-slices. Because of this, we use only the Sp and Wp results to 
extract scattering amplitudes. The data show that the interaction in the D + C channel 
is attractive, consistent with the result derived in the Appendix. 

Our results for scattering amplitudes are presented in Table 3. We illustrate our 
procedure for obtaining these results using T(S) as an example. For each jackknife sample 
we first average the fit value for 6E(S) from the Sp and Wp sources, and then solve the 
cubic polynomial given in Eq.(2.1) for T(S). The central value and the error are given by 
the jackknife procedure, regarding the 35 data points as statistically independent. Within 
the same jackknife procedure we also calculate T(A), and extract T q = {T(S) — T(A))/2 
and T g = (T(S) + T(A))/2. 

When solving Eq. (2.1) for T(S) or T(A) we monitor the effect of the l/L 4 and 1/L 5 
terms (using the values of given in Table 1). These turn out to be, respectively, ~ 31% 
and 8% of the leading term. This suggests that the error introduced by truncating the 
series is a few percent in the scattering amplitudes. This is smaller than the statistical 
errors, which are approximately 10%. 

To test the current algebra predictions, we calculate the combinations 4/^T. These 
are included in Table 3, and shown in Fig. 8. The errors are calculated assuming that 
/tt and T are uncorrelated; we have checked on subsamples that the errors are similar if 
correlations are included. In the chiral limit the following relations should hold: 4/^T(S") = 
1, T(S) = T q = —T(A) and T g = 0. Our results show that, within our errors, the second 
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relation holds even at the relatively large quark masses we have used, and Af^T lies between 
0.76 — 0.88. There is a small increase in T between n = 0.154 and 0.155, but because of the 
size of the statistical errors we cannot conclude if this is related to the \jm^ divergence 
expected for Wilson fermions in the chiral limit. Also, the errors are too large to extract a 
value for the gluon exchange amplitude T g . All we can say is that it is much smaller than 
the quark exchange amplitude. 

It is interesting to compare our results with those obtained using staggered fermions 
on the same lattices [4]. For technical reasons, we were able to calculate only the quantity 
AQ = T q — 2T q T g (cim 7T / AttL)) with staggered fermions, and not T q and T g separately. 
If we use the values for T g obtained here, however, then we find that AQ = T q to good 
approximation. To compare the results, we note that the pion mass and f K match for 
the following parameter values: staggered m q = (0.02 + 0.03) with Wilson n = 0.154, and 
staggered m q = (0.01 + 0.02) with n = 0.155 [11]. The staggered results at these two 
masses are T q AQ = 67(8) and 78(6) respectively, which agree within errors with the 
results obtained here. It is reassuring to find that the two formulations, each with their 
separate technical problems, yield mutually consistent results. 



5. Expected behavior of lZ(t) with Wilson Fermions 

As shown in Table 2 and Figs. 4-7, different source and sink operators give consider- 
ably different values for Za,s- This variation is more marked for smaller quark mass. It 
turns out that the values for Za,s can be understood semi-quantitatively using a simple 
model, as we explain in this section. This model is similar to that used in the analysis of 
our staggered fermion data [4]. 

We begin by imagining that the source creates two pions each having p = 0. The 
pion operators could be local or smeared, but should have a finite extent that is much 
smaller than the lattice size. Also, the same source operators are used in the numerator 
and denominator of Eq. (3.4). In this case we expect that Za,s 1 as X -> oo, because 
the lightest two pion state differs from two independent, zero momentum pions by terms 
which vanish as L — > oo. Thus the Z-factors in the numerator and denominator of Eq. 
(3.4) cancel. Assuming Za,s = 1? an d using Eqs. (3.5) and (3.6), the behavior of the 
direct and crossed contractions is 

D(t) = l-^±^ + 0(t 2 ), 

F F ^ 
C(t) = _^^ + 0(t 2 ), 



for t large enough that the contributions of excited states have died away. What is impor- 
tant here is that only D(t) is non-zero when extrapolated back to t = 0; the constant term 
corresponds to the two pions propagating from source to sink without interactions, which 
is only possible in the Direct channel. In the Crossed channel, at least one quark exchange 
interaction is required. This gives rise to the term linear in t, since the interaction can 
occur at any time. The term linear in t in D(t) is due to the gluon exchange interaction. 

In practice Za and Zs differ from unity, for two reasons. First, the two pion state is 
altered by interactions. This gives corrections proportional to 1/L 2 [4], which we assume 
are small and ignore. Second, our sources are not two independent pion operators each 
having p = 0, but rather a single wall or Wuppertal quark source. This gives large 
corrections to Za,s, due to the overlap, which do not vanish as L — > oo, and it is these 
which we estimate. 

In our set up a state of two quarks and two antiquarks, all in close proximity to 
one another, is created rather than two separate pions. We attempt to pair the quarks 
and antiquarks into pions by making two color singlets, each with pseudoscalar quantum 
numbers. But, by performing a combined color and Dirac Fierz transformation, we find 
that we are also creating, with non- vanishing amplitude, two pions with the opposite qq 
pairings. Explicitly, the Fierz transformations are 



where V ®V = (V'iTsV^XV'sTsV^) and A ® A = (V'i7475V'2)(V'37475V'4), parentheses 
implying spin and color traces. The sign due to fermion exchange is not included in these 
Fierz identities, since we do not include the sign in our definitions of D(t) and C(t). The 
Fierz relations hold for both wall and Wuppertal sources, in the latter case because all 
the products of link matrices begin at the same site. We have shown only the V ®V and 
A <E> A parts of the Fierzed combinations because the correlators of these operators have 
the dominant contribution. Other tensor structures give no contribution for two separated 
pion sources in infinite volume, and thus give contributions here that are suppressed by 
powers of The same is true for operators consisting of two color octets. The magnitude 
of these neglected terms can be significant, especially for Wuppertal sources, as shown by 
the difference between our data and the estimates presented below. 




(5.2) 
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There is no large Fierz contribution from the sinks, since these do consist of two 
independent pion operators, each having p = 0. The only contribution occurs when the 
two operators overlap, and is suppressed by powers of 1/L. 

The Fierz relations mean that our Crossed contractions contain a part in which the 
quarks have already been exchanged before we can identify the state as one with two pions, 
so that no subsequent quark-exchange interaction is necessary. This leads to a constant 
term in the Crossed contraction. The Fierz contributions to the Direct contraction do not, 
however, affect the constant term, for there must be an additional quark exchange inter- 
action to bring the quarks and antiquarks back to their original pairings. This discussion 
motivates the following assumptions for the constant terms 

D(t = 0) p ,a ~ 1 , 
C(t = 0) P « ^(1 + C 2 A /C 2 P ) , (5.3) 
C(t = 0) A « ^(C P /C A + 1) . 

The subscript indicates the type of correlator, and the constants C a-, Cp are the amplitudes 
for creating single pions with p = using the operators A4 and V respectively. The ratio 
Cp/Ca is 2.5 and 3.1 for Wuppertal sources, and 1.6 and 1.9 for wall sources, at n = 0.154 
and 0.155, respectively. Using these values, we can calculate the Z's using 

Z s « D(t = 0) - C(t = 0) , Z A « D(t = 0) + C(t = 0) . (5.4) 

The predictions are collected in Table 4. They give a good semi-quantitative represen- 
tation of the data for Za,s m Table 2. In particular, we can understand the small value of 
Zs for the Sa operators as being due to a large cancellation between the Direct contraction 
and the Fierz contribution to the Crossed contraction, the latter being enhanced by the 
large ratio Cp/Ca for Wuppertal sources. This cancellation is most likely why the signal 
is so noisy in this channel. 

6. Conclusions 

We find that it is straightforward to calculate the finite volume energy shift for chan- 
nels not involving qq annihilation. The calculation is much less involved using Wilson 
fermions than that we carried out with staggered fermions [4]. We are able to work on 
a lattice of modest size (L 1.6fm) because the interactions in the channels we consider 
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are relatively weak. From the energy shifts we extract the quark exchange amplitude, 
and place a bound on the gluon exchange amplitude. Our results are consistent with the 
predictions of current algebra, on the other hand the quarks used in the calculation are 
not light enough to expose the expected artifacts due to the breaking of chiral symmetry 
by Wilson fermions. 

It is important to extend this work to smaller quark masses, where the divergence 
in T due to chiral symmetry breaking should show up. Furthermore, the result should 
be checked on a larger volume to verify that the asymptotic form of the finite volume 
dependence can be used. 

Acknowledgments 

The 16 3 x 40 lattices were generated at NERSC at Livermore using a DOE allocation. 
The calculation of quark propagators and the analysis has been done at the Pittsburgh 
Supercomputing Center, San Diego Supercomputer Center, NERSC and Los Alamos Na- 
tional Laboratory. We are very grateful to Jeff Mandula, Norm Morse, Ralph Roskies, 
Charlie Slocomb and Andy White for their support of this project. AP thanks Los Alamos 
National Laboratory for hospitality during the course of this work. SRS is supported in 
part by DOE contract DE-AC05-84ER40150 and by an Alfred P. Sloan Fellowship. 

Appendix A. Correlation Inequalities for Scattering Amplitudes 

One can derive rigorous inequalities among correlation functions for vector-like gauge 
theories such as QCD. The basis of such inequalities is the positivity of the measure in the 
Euclidean path integral: 

[dA^][dtp][dtp] exp[-S g aU g e - Sf ermion ] > . (A.l) 

This property has been exploited, both on the lattice [12] and in the continuum [13], to 
derive inequalities among 2— point correlation functions. As a result it was shown that 
the pion is the lightest meson. Here we apply the same arguments to 4— point correlation 
functions in finite volume to derive constraints on tt — tt scattering amplitudes. The only 
other work extending the derivation of inequalities to higher order correlation functions 
that we are aware of is Ref. [14], where it is shown that the pion wavefunction is largest 
at r = using a 4— point correlation function. 
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We consider the 4— point correlation functions corresponding to the tt — tt scattering 
at threshold. Let the sources for the two pion be at time t = and the two pion sinks 
be at time t = T . We take all four pion operators to be point-like with Dirac structure 
V'TsV' an d zero 3— momentum. In terms of the quark propagator G(x } 0; y, T), the Direct 
and Crossed correlators are: 

D(0,T;p = 0) = ]T (Tr(G(f,0;y,T)Gt(f,0;y^ , 

x ,y, z,w 

C(0,T;p = 0) = ]T (Tr(G(f,0;y,r)Gt(£,0;y,r)G(£,0;^,r)Gt(f,0;^,T))) , 

x ,y, z,w 

(A.2) 

where the trace is taken over the color and spin indices. 
Using the Schwarz inequality 

(fP) > \(f)\ 2 , (A.3) 

we get the relation 

D(0,T;p = 0) > [P(0,T;p = 0)] 2 (A.4) 
where the zero 3— momentum pion correlator, 

P(0,T;p = 0) = J](Tr(G(f,0;y,T)Gt(f,0;y,T))) , (A.5) 

is by itself positive definite. This inequality implies that the contribution of this diagram 
to the two pion interaction is attractive. 

These results can be generalized to other Dirac structures. In fact the interaction 
between any two mesons, e.g. two rhos, is attractive in the Direct channel. Note that 
the derivation of this inequality did not depend on the volume of the system. This is 
analogous to the zero temperature mass inequalities of Ref. [12] [13] which can be used at 
finite temperature to give relations between hadronic screening lengths. 

The crossed correlator given in Eq. (A.2) can be rewritten as 

C(0,T;p = 0) = ]T(Tr([]TG(f,0;y,T)G^^ , 

x ,z y w 

(A.6) 

and is therefore also positive. Combining this fact with the inequality in Eq. (A.4) shows 
that the scattering amplitude in the flavor antisymmetric channel, corresponding to the 
combination D + C, is 

D(0,T;p = 0) + C(0,T;p = 0) > [P(0, T;p = 0)] 2 . (A.7) 
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This inequality on the correlators implies that the exponential fall-off with time in this 
channel is slower than that for two non-interacting pions. It follows that the interaction 
energy SE in the antisymmetric channel is negative, i.e. the scattering length is positive. 
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p 




Lattice 


N con f 


m w 


U 


6.0 


0.154 


16 3 x 80 


35 


0.364(6) 


0.057(3) 


6.0 


0.155 


16 3 x 80 


35 


0.297(9) 


0.055(3) 



Table 1. Summary of results from Ref. [8] needed in this calculation. f K is normalized 
such that the experimental value is f K = 93 MeV , and is obtained using the mean field 
improved value 0.77 for the axial current renormalization constant. 





D + C 


D-C 


Correlator 


Fit 


Z A 


-E A 


Fit 


Z s 


Es 


k = 0.154 


s P 


9-16 


1.19(4) 


0.020(2) 


8-25 


0.93(7) 


0.022(6) 


W P 


8-16 


1.10(2) 


0.022(3) 


8-15 


0.92(2) 


0.018(4) 


Sa 


8-16 


1.77(10) 


0.019(7) 


8-20 


0.46(7) 


0.019(10) 


w A 


8-15 


1.30(5) 


0.025(6) 


8-26 


0.78(7) 


0.016(9) 


k = 0.155 


s P 


8-17 


1.23(6) 


0.022(6) 


8-18 


0.87(6) 


0.022(8) 


W P 


8-18 


1.14(3) 


0.027(3) 


8-15 


0.92(3) 


0.024(5) 


Sa 


8-15 


2.07(17) 


0.020(9) 


8-20 


0.27(11) 


0.023(25) 


w A 


8-15 


1.49(9) 


0.023(5) 


7-20 


0.67(8) 


0.017(13) 



Table 2. Results for the amplitude and energy shifts obtained from fits to lZ{t) of cor- 
relators for the D ± C channels. The four kinds of correlators, Sp } Wp } Sa and W a are 
described in the text. We also give the range of time-slices over which the fit is made. 
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k = 0.154 


k = 0.155 


T(S) 


58.7(8.7) 


69.6(10.7) 


T(A) 


-61.8(6.5) 


-73. 0( 6.6) 


T q 


60.2(4.6) 


71. 3( 6.6) 


T g 


-1.6(5.0) 


-1.7( 6.0) 


±flT(S) 


0.76(14) 


0.84(16) 


AflT(A) 


-0.80(12) 


-0.88(13) 




0.78(10) 


0.86(12) 




-0.02( 7) 


-0.02( 7) 



Table 3. Final results for the scattering amplitudes. 



Correlator 


Z a 


Zs 


k = 0.154 


s P 


1.10 


0.90 


W P 


1.12 


0.88 


Sa 


1.60 


0.40 


W a 


1.30 


0.70 


k = 0.155 


Sp 


1.09 


0.91 


W P 


1.11 


0.89 


Sa 


1.88 


0.12 


W a 


1.38 


0.62 



Table 4. Model predictions for the intercepts of lZ{t). 
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(a) 




Fig. 1. The four different contractions that contribute to the two pion correlator: a) 
Direct, or gluon exchange, channel (-D), b) Crossed, or quark exchange, channel (C), c) 
Single Annihilation and d) Double Annihilation. The diagrams also correspond to the 
amplitudes that contribute to tttt scattering. 
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Fig. 2 The effective energy shift 6E(t) using Sp correlators at n = 0.154, for both the 
1 = 2 (S) representation (D — C), and the A representation (D + C). 
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Fig. 3 As in Fig. 2 but for the Wp correlators. 
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4. lZ{t) using Sp correlators at n = 0.154. 
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6. lZ{t) using Sa correlators at n = 0.154. 
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Fig. 8. plotted versus m^/m 2 p to test the chiral behavior. The data for T q has been 

displaced by —0.05 along the x-axis for clarity. 
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